Faraday waves in elongated superfluid fermionic clouds 



OO 
O 

o 

(N 
Os 



5-H 



I 

C 

O 

o 



> 
Os 

in 

06 
o 

OO 

o 



X 



P. Capuzzi^B and P. Vignolo 2 -Q 

1 Consejo National de Investigaciones Cientificas y Tecnicas and Departamento de Fisica, 
FCEyN, Umversidad de Buenos Aires, 1428 Buenos Aires, Argentina 
2 Institut Non Lineaire de Nice, Universite de Nice-Sophia Antipolis, 
CNRS, 1361 route des Lucioles, 06560 Valbonne, France 

We use hydrodynamic equations to study the formation of Faraday waves in a superfluid Fermi 
gas at zero temperature confined in a strongly elongated cigar-shaped trap. First, we treat the role 
of the radial density profile in the limit of an infinite cylindrical geometry and analytically evaluate 
the wavelength of the Faraday pattern. The effect of the axial confinement is fully taken into account 
in the numerical solution of hydrodynamic equations and shows that the infinite cylinder geometry 
provides a very good description of the phenomena. 
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I. INTRODUCTION 

The term Faraday waves refers to regular surface 
fringes which are excited by a vertical oscillatory mo- 
tion in a nonlinear liquid The time-dependent spa- 
tially uniform driving leads via the nonlinear interaction 
between the wave excitations of the system to an insta- 
bility towards the formation of spatial structure. Since 
the Faraday discovery, the phenomena of pattern for- 
mation and, more generally, of parametric amplification 
have been studied in diversified contexts, including con- 
vective fluids, nematic liquid crystals, nonlinear optics 
and biology Q- Recently, Faraday waves have been ex- 
perimentally created in a trapped atomic Bose-Einstein 
condensate (BEC) by a periodic modulation of the trans- 
verse confinement Nonlinearity in ultracold gases can 
be driven either by varying the scattering length, e.g., 
via Feshbach resonances, as proposed in this context by 
Staliunas et al. [4] , or by varying the trap parameters 
as experimentally realized by Engels et al. [3j]. Several 
phenomena connected to parametric amplification of ex- 
citations has been theoretically investigated in atomic 
BECs [1, H, 0] and very recently, the properties of Fara- 
day waves in confined BECs has been also explored [1, Q . 
In all these cases, the periodic modulation of the nonlin- 
earity leads to a parametric excitation of sound waves in 
the transverse direction respect to the modulation. 

In this paper, we study the Faraday pattern formation 
in a superfluid Fermi gas in the BCS-BEC crossover at 
zero temperature. The physics of Faraday waves in a 
superfluid Fermi gas is richer than in an atomic Bose- 
Einstein condensate because the nonlinearity may also 
drive the microscopic features of the superfluid as de- 
scribed by its equation of state. This leads to strong 
variations in the pattern along the crossover. 

We adopt the experimental conditions of Engels et al. 
where the radial confinement of an elongated trap is 



periodically modulated in time. By using a hydrody- 
namic description of the superfluid Fermi gas, we first 
treat the case of a infinite cylindrical gas to obtain an 
analytical expression for the spatial modulation of the 
Faraday pattern as a function of the driving frequency VL, 
for small frequencies Q <C u>±, luj_ being the transverse 
trap frequency. Then we numerically solve the three- 
dimensional hydrodynamic equations taking into account 
the axial confinement and investigate the formation of 
Faraday waves at higher frequencies (Q ~ 2lo±_) where 
the fringes in the density profile can be experimentally 
visualized in cigar-shaped traps. Finally we compare the 
analytical prediction for the infinite cylinder with the nu- 
merical results and find that our analytical formula pro- 
vides a good estimate of the pattern modulation even at 

n ~ 2u±. 

The manuscript is organized as follows. Section |TT] in- 
troduces the hydrodynamic description of the gas and 
derives the low frequency solution of the spatial pattern 
formation. In Sec. Mil we present the numerical results 
for the dynamics in a cigar-shaped trap and compare 
them with our analytical prediction. Section IIVI offers a 
summary and concluding remarks. 



II. FARADAY EXCITATIONS IN A CYLINDER: 
HYDRODYNAMIC DESCRIPTION 

Let us consider a superfluid Fermi gas at T — con- 



fined by a cylindrical potential Vj_(r) 



' r /2 with 



r = (x + y ) and m the atom mass. The dynamics of 
the density profile n(r,t) and the velocity field v(r, f) is 
described by superfluid hydrodynamic equations, namely, 
the continuity equation 
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where fi[n) is the equation of state (EOS) of the gas 
which contains the microscopic details of the superfluid. 
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A solution of the dynamics for a time-dependent u>±_ 
can be obtained by relying on a scaling Ansatz. Follow- 
ing Ref. [13] we introduce the scaling parameter bit) and 
write 



n(r,z,t) = 



1 



b(ty 



:no(r/b(t)) 



and 



(3) 



where no(r) is the stationary solution of the Euler equa- 
tion at t = 0, i.e., the Thomas- Fermi (TF) profile defined 
by n(no(r)) = [p, — V±(r)] with /2 the chemical potential. 
The scaling Ansatz in Eq. ((3]) satisfies the continuity 
equation for any form of bit) and is an exact solution 
of the Euler equation if the EOS follows a power law, 
fi(n) = Cn 1 . In this case the scaling parameter obeys 
the differential equation 



b + LU 2 ± it)b- 



= 0. 
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Notice that if the EOS is not a power law, one can still 
introduce an effective exponent 7 and exploit the same 
procedure [l2| to find an approximate solution of the dy- 
namics. 

To excite Faraday waves, we modulate the transverse 
trap frequency as 



u±(t) 



.(1 + e sin fit) 
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where Q is the driving frequency and e is a small pa- 
rameter, e -C 1. For fl <C Eq. fU has the solution 

El 



b(t) ~ 1 - a sin(flt) 



(6) 



with a — e/(-f + 1 — fl/(2uj±)) ~ e/fa + 1) in agreement 
with the expression found in Ref. [8| for an ultracold 
Bose gas (7 = 1). For larger fi, the solution of Eq. ((4]) 
represents a forced breathing mode which does not adia- 
batically follow the forcing frequency. For instance, when 
the driving frequency is close to the natural breathing 
mode ft — y/2~f + 2ui±, the system exhibits the resonant 
solution 



bit) ~ i + EUJ1 - t cos(nt) 



(7) 



valid for small t. For larger t, the resonant solution can 
be seen to exponentially increase. 



A. Excited states 

The modulation of the nonlinearity excites sound-like 
phonons in the z-direction which, in turn, lead to the in- 
stability of the scaling solution J3]) . These excited states 
breaking the axial symmetry can be investigated by a 
linear analysis, i.e., by setting n(r,z,t) = n$(r jb)jb 2 + 
<5n(r, z, t) and linearizing Eqs. |lj and (|2|) in 5n we arrive 
the following equation 
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that is valid for small velocity fields, i.e., SI <C Equa- 
tion |5]) is a sound-wave equation for the instantaneous 
density. Using the solution Eq. ([6|) and developing to 
first order in a, Eq. ((SJ) gives rise to 
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where R± = y/2p,/muj 2 _ is the TF radius of the gas. 
The right-most term in Eq. ([9]) proportional to sin(Slt) 
behaves like a source for excitations in the linear re- 
sponse regime. Expressing Sn in its Fourier compo- 
nents, Sn — J2j Snj(r, z) cos(ujjt + tfij), one can easily 
demonstrate that the lowest-frequency Fourier compo- 
nent that is resonant with the perturbation corresponds 
to wo = fi/2. This is a general result for dynamical 
systems governed by Mathieu-type equations [13j as Eq. 
© . We keep only the resonant term in Eq. §§§ and make 
the Ansatz (5no(r, z) = Sno(r) cos(gz), corresponding to 
counter propagating axial phonons with momenta ±q. To 
find an estimate of q, we assume that fl is small enough 
to consider only the sound mode Sno and therefore take 
Sno{r) = (dfx/dn) 1 evaluated at a = 0, i.e., the sound 
limit solution for q — > and ft — > at t = (see Ref. 
[Til]). By integrating in the radial plane one obtains 
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which, for a power-law EOS leads to 



q 




(11) 



(12) 



neglecting second-order terms in a. The above expression 
can be directly linked to the sound velocity c s [14] in an 
elongated superfluid by q = Q/(2c s ). 

By introducing the wavelength of the Faraday pattern 
d = 2n/q it is possible to evaluate the number of fringes 
M in an elongated superfluid confined by an harmonic 
potential V ex tir, z) — V±(r) + muj 2 z 2 /2. Given that the 
axial extent of the density can be written as L z = 2R±/X 
where A = uj z /uj± and neglecting edge effects, Eq. (TT^)) 
yields 



d 7rA V 27 dJj_ 



(13) 



For 7 = 1, M is in good agreement with the experi- 
mental data of Engels and coworkers [H even for large 
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s and il, where the procedure we used to derive it is 
not fully justified. Equation (| 1 3[) predicts that the num- 
ber of Faraday fringes varies by a factor of 2/y5 in the 
BCS-BEC crossover and thus it could be easily observed 
in current experimental setups with strongly elongated 
traps HI, [l6| , as discussed in the following section. 



III. NUMERICAL RESULTS 

The numerical solution of the superfluid hydrodynamic 
equations allows to fully take into account the role of the 
axial confinement and to get over the approximation of 
small frequencies. To compute the dynamics we first map 
the superfluid density n(r,t) and velocity field v(r, t) to 
the complex field ip(r,t) = y/ne lS , with v = hX7S/m. 
In this way the dynamical equations can be cast in a 
Gross-Pitaevskii-like equation of motion for ip. We then 
discretize this equation in the (r, z) plane and solve it 
using a 5th-order Runge-Kutta method [l7|. A mesh of 
200 x 400 points in the (r, z) plane and a time step dt = 
10~ 4 /w_l are usually enough to compute the evolution 
until t = 500/w^. 

We consider a Fermi gas with N = 5 x 10 5 40 K 
atoms trapped in a cylindrically symmetric trap with 
lu± = 160.5 Hz and uj z — 7 Hz, and fix the modula- 
tion parameter e = 0.05. The superfluid is described by 
the EOS developed by Hu et al. [3 and parametrized 
in terms of l/(fcp-a) where Uf is the Fermi wavevector 
kp = (37r 2 n) 1 / 3 and a is the s-wave scattering length. 
If not differently specified, in the following we refer to a 
continuous modulation of the transverse trap frequency 
at Q = 2luj_. 

From the numerical solution of the density profile 
n(r, z,i) we have computed the axial and radial RMS 
widths, and the column densities nuy(z,t) = J dyn(x = 
0, y, z, t). These are shown in Figs. [JJ3] together with 
their Fourier transforms for the case of the Fermi gas 
in the BCS regime (a/ao = — 10 _1 , ao being the Bohr 
radius). The dynamics of the width Ar (Fig. [JJ is dic- 
tated by the driving frequency coupled to the transverse 
breathing mode, which for the gas in the BCS limit is 
to = ^/TO/3 This reflects the competition between 
the forcing and the natural breathing mode. Further- 
more, within the scaling approach of Sec. HU the axial 
width evolves in time with bit), i.e., Ar(t) = Ar(0) b(t). 
Indeed, the numerical solution of the ordinary differen- 
tial equation ([4| reproduces very well the behavior of Ar 
for the axially confined gas before the Faraday instability 
sets in. As noted earlier, we would like to emphasize that 
due to the vicinity of the natural breathing mode and the 
driving frequency, Ar does not simply adiabatically fol- 
lows lo± as clearly shown in the figure. 

The time evolution of the axial width Az and its 
Fourier transform are instead illustrated in Fig. [2j Sim- 
ilarly to the radial width, the evolution of Az is domi- 
nated by the forcing frequency £1 and the axial breathing 
mode, that is y/12/5uj± in the BCS case. The excita- 
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FIG. 1: Time evolution of the transverse width Ar for a Fermi 
superfluid in the BCS limit (see text). The top panel depicts 
Ar (in units of aho = \J Ti / (muj ±)) as a function of t (in units 
of loJ 1 ) and the bottom panel depicts its Fourier transform 
(in arbitrary units and logscale) as function of u/ui_l. 
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FIG. 2: Idem Fig. \T\ for the axial width Az. The Fourier 
transform in the bottom panel is displayed in arbitrary units 
and logscale, and the frequency u> in units of uo z . 



tion of the axial breathing mode is a consequence of the 
coupling of the radial and axial directions imposed by 
the three-dimensional trap and should not be expected 
in the infinitely long cylinder discussed in Sec. [Til 

Although Ar and Az show no evidence of the fre- 
quency f2/2 in their spectra, the column density at the 
center of the trap (Fig. [3]) clearly exhibits it. The relative 
importance of this frequency in the spectrum of ?iid is 
accentuated with time as a result of the formation of the 
Faraday wave, this effect can be seen by analyzing the 
Fourier transform for different time intervals. 
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FIG. 3: Time evolution of the column density at the trap 
center n? D = nm(z — 0). The top panel shows n? D (in units 
of ) as a function of time t (in units of u>±). The bottom 
panel displays its temporal Fourier transform (in arbitrary 
units and logscale) as function of uj/uj±. 

A. Onset of Faraday wave formation 

The Faraday instability becomes visible when the mo- 
mentum distribution develops components at the mo- 
menta ±q corresponding to pairs of counter-propagating 
axial excitations of frequency ft/2. In Fig. Q]we display 
the formation of a Faraday wave in the density profile of 
a gas in the weakly BCS limit with a = — 10 _1 ao. Af- 
ter approximately 50 driving periods, the appearance of 
transverse fringes along the axis of the cigar is clear and 
we can say that the Faraday instability has been fully set 
in. For longer times, as the Faraday waves increase their 
amplitude, a large number of modes interacts and the 
Faraday pattern develops several complicated spatial and 
temporal structures, resembling those of spatio-temporal 
chaos as expected in confined hydrodynamic flows Q. 
This regime cannot be accounted for by the simplified 
analysis of Sec. |TT] and requires a more demanding nu- 
merical effort. 

The onset time of the Faraday waves varies through 
the BCS-BEC crossover. The growth rate of the spatial 
pattern can be evaluated by computing the fraction 

^'^{jM^k) (14) 

where ho(z,t) is the scaled density extracted from the 
numerical solution by fitting a TF profile n(z) — A{jl — 
rrvJ^z 1 /2) v to the integrated density nxD(z,t) at each 
time t, and Sn = riio{z,t) — ho(z,t). The magnitude Sn 
represents the normalized deviation of the density respect 
to the scaling solution and can therefore be related to 
the amplitude of the Faraday wave. The evolution of 
<5jv is illustrated in Fig. for three values of a in the 
BCS, unitarity and BEC regimes. In the BCS and in 
the unitarity limits the growth rates are comparable as 



FIG. 4: (color online) Density plots n(x, y = 0, z, t) of a Fermi 
gas in the BCS regime at several times t during the transverse 
frequency modulation. 




co ± t 

FIG. 5: Fraction Sn (in logscale) as a function of the time t 
(in units of wj ) for a/ao = — 10 _1 (circles) , — 10 5 (triangles), 
and 10 3 (squares). The empty symbols corresponds to the 
continuous driving while the full ones to the 3-cycle driving. 
The lines are guides to eyes for better visualizing the growth 
rates. 

the excitation spectra in both regimes are the same. In 
the BEC limit, instead, the onset is more rapid since the 
driving frequency ft = 2lu± corresponds to a collective 
mode of the cloud. 

Faraday waves can be also created by driving the trans- 
verse confinement for only few cycles [3| . However, in this 
case, the growth rate of the pattern is modified. For in- 
stance, we have found that for a 3-cycle perturbation, 
the rate for a BCS and a unitarity superfluid Fermi gas 
increases approximately by a factor of 2-3, while in the 
BEC limit it decreases by a factor of 1.5 approximately. 
For our choice of ft, the driving frequency and the breath- 
ing modes are in competition in the BCS and unitarity 
limits and the continuous excitation may delay the for- 
mation of the pattern. In the BEC limit, ft = 2u± corre- 
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sponds to the breathing mode and the onset of the insta- 
bility is more rapid by continuously exciting the system 
at resonance. Note that the total growth time of a Fara- 
day wave may also depend on the numerical noise caused 
by the discrete space-time grids, but the growth rates in 
the linear regime should not be modified. 



B. Faraday wavevector through the crossover 

A further observable which depends on the crossover 
is the wavevector q of the Faraday pattern, as already 
outlined in Sec. [Ill Experimentally, q can be deduced by 
measuring the interfringe distance d at the center of the 
cloud. The interfringes obtained in the numerical sim- 
ulation can be observed in Fig. [6] where the integrated 
density ti^d has been plotted together with the topview 
density plots for a/ao = — 10 -1 ,— 10 5 , and 120 corre- 
sponding to the BCS, unitarity and BEC regimes, re- 
spectively. It is worth stressing that at a = 120ao, the 
transverse size displays an impact-oscillator behavior as 
observed for an atomic BEC at the breathing frequency 
[H, Q . The resulting strong modulation of the nonlinear- 
ity speeds up formation of the Faraday waves and impose 
severe numerical limitations. 

Through the crossover, from the BCS to the BEC lim- 
its, at fixed number of fermions, both the size of the cloud 
and the interfringe d decrease since they are both propor- 
tional to yfp.\ however, d decreases less rapidly because of 
the factor y/j/Jl + 7) with a corresponding slight modi- 
fication of the fringe number. As a result in the BCS and 
unitary regimes the number of fringes are the same as 
can be seen in Fig. 6, i.e., 16 for a = — O.lao and — 10 5 ao, 
while in the BEC regime (a = 120ao), the number of 
fringes diminishes by a factor of 2/y/E, corresponding to 
14 fringes approximately. In more elongated traps as that 
used in the ENS experiments [l6j with A ~ 0.18, the vari- 
ation of Af through the crossover can be more important, 
from M = 41 in the BCS limit to M = 37 in the BEC 
limit at f2 = 2lu±. 

In Fig. [7] we compare the Faraday wavevector q — 2-k / d 
as observed in the numerical simulation with that pre- 
dicted for a soundlike excitation at f2/wj_ = 2 and 1.5. 
Mimicking a real experiment, d has been measured at the 
center of the density profile and the error bars have been 
evaluated from its variance during several time steps. 
Equation pT|) instead has been evaluated by using the 
EOS deduced by Hu and coworkers [HI]. The results show 
that Eq. (fTTj) provides a good estimate of the function 
q(£l) through all the crossover in particular for Q = 2lu±. 
On the other hand, for Q = 1.5u>± the numerical result 
for a — 10 3 ao lies well above the prediction. Indeed, by 
analyzing the FFT of n® D in this case, we found that the 
transverse modulation at Q = 1.5cj_l also excites the col- 
lective mode of the BEC at f2 = 2oj± similarly to what 
was found experimentally in Ref. j|. As a consequence, 
the mode at f2 = 2lu±_ dominates the generation of longi- 
tudinal excitations even when modulating at f2 = 1.5uj± 




FIG. 6: (color online) Density plots and integrated density 
tjid as function of position (in units of dho) after the Faraday 
instability has been set in. Top, middle, and bottom panels 
correspond to a/ao = 120, —10 s , and — 10 -1 , respectively. 



and thus the numerical result should be compared with 
the prediction at U, = 2luj_, as shown by a cross in Fig. 
[7] For comparison, we have also evaluated the Faraday 
wavevector predicted by Eq. (11) using the mean-field 
BCS EOS, /i(n) = h 2 k 2 F /2m(l + ^k F a). It is worth- 
while remembering that this EOS predicts the collapse of 
the fermionic mixture for strong attractive interactions, 
as indicated by the divergence of q shown in dashed line 
in Fig. [3 



IV. SUMMARY AND CONCLUDING 
REMARKS 



In summary we have studied the formation of Fara- 
day patterns in an elongated supcrfluid Fermi gas in the 
BCS-BEC crossover. We simulated an experiment where 
nonlinearity is driven by modulating the transverse con- 
finement and compared the numerical results with an an- 
alytical estimate derived neglecting the axial confinement 
and for low frequencies. The peculiarity of a superfluid 
Fermi gas respect to an atomic Bose-Einstein conden- 
sate is that the nonlinear term controls the microscopic 
features of the superfluidity and that the pattern mod- 
ulation depends on the crossover itself. Moreover the 
pattern growth rate varies through the crossover as a 
consequence of the different spectra. Our numerical and 
analytical predictions may be observed in current exper- 
iments with superfluid Fermi gas in strongly elongated 
traps [HGit. 
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FIG. 7: Comparison of the ratio q/Q between the observed 
Faraday wavevector q (in units of aho) and the driving fre- 
quency Q (in units of lu±) with that predicted for a soundlike 
excitation (Eq. (JTTJ) ) as functions of l/(kpa), where k° F is 
the Fermi momentum of the noninteracting gas at the trap 
center. The error bars have been estimated from the vari- 
ance of d during several time steps. The cross corresponds 
to a result obtained by modulating at Q/ui± = 1.5 corrected 
by taking into account the induced excitation at n = 2lu±_ 
(see text). The dashed line corresponds to the prediction of 
Eq.JTTJ with the BCS mean-field EOS. 
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